arXiv:1507.02885vl [cond-mat.mes-hall] 10Jul2015 
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We demonstrate theoretically that charge transport across a Josephson junction, volt age-biased 
through a resistive environment, produces antibunched photons. We develop a continuous-mode 
description of the emitted radiation field in a semi-infinite transmission line terminated by the 
Josephson junction. Within a perturbative treatment in powers of the tunneling coupling across the 
Josephson junction, we capture effects originating in charging dynamics of consecutively tunneling 
Cooper pairs. We find that within a feasible experimental setup the Coulomb blockade provided by 
high zero-frequency impedance can be used to create antibunched photons at a very high rate and 
in a very versatile frequency window ranging from a few GHz to a THz. 
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Introduction Photons from ordinary thermal sources 
have a tendency to bunch together, and the first con¬ 
trolled generation of single photons was only performed 
in 1974 [l|. Since then, single photons have been used 
to explore fundamental aspects of quantum mechanics, 
such as the interference fringes of single particles 0 and 
Bell’s inequalities A good single-photon source can 
be characterized by the vanishing probability of detecting 
two photons at the same time from its output. This prop¬ 
erty is called anti-bunching Q , and is also a sign that the 
field is non-classical In addition to being of funda¬ 
mental interest, the capability to create single propagat¬ 
ing photons is an indispensable tool for many quantum 
information applications, including quantum key distri¬ 
bution [^. 

In the optical domain, the workhorse for many single¬ 
photon experiments has been parametric down conver¬ 
sion [^. Here, photon pairs are generated at random 
times, and one of the two photons can be used to herald 
the other one . Other type of sources for photons in the 
visible spectrum are based on photoluminescence in sys¬ 
tems such as single ions, single molecules, semiconductor 
quantum dots and diamond color centers [11^. 

In the microwave domain, superconducting quantum 
circuits have been used to build single photon 

sources [la - H^ . In most of these realizations, a two-level 
system is first excited and consequently spontaneously 
emits a microwave photon into a transmission line. The 
emission rate can often be controlled in-situ, allowing for 
waveform shaping plrl - fT^ . The use of the field reflected 
from a single artificial atom, which is indeed perfectly 
anti-bunched, has equally been demonstrated [20l - l^ . In 
addition to fundamental quantum optics experiments, 
single photon sources in the microwave domain might be 
useful for metrological purposes due to their well defined 
output power and, in combination with single photon de¬ 
tectors, for quantum nondemolition measurements due to 
their ideal amplitude squeezing. 

In this article, we study theoretically an alternative 


type of single-photon source, based on the Coulomb 
blockade of charge transport across a Josephson junction 
in series with an electromagnetic environment (23l - [^ . 
Our proposed setup has a very simple operating prin¬ 
ciple and promises to be a bright, robust and versatile 
’photon gun’. 

Anti-bunched photons are created from the ap¬ 
plied voltage, making use of inelastic (photon-assisted) 
Cooper-pair tunnelling [2^ - 1^ and the long charging 
time of the junction. When biased below the supercon¬ 
ducting gap, the applied voltage defines the frequency 
spectrum of the emitted photons through the Josephson 
frequency, a;j/27r = 2eV/h. This frequency is not limited 
to the plasma frequency of the Josephson junction and 
the energy of the emitted photons is, in principle, only 
limited by the gap of the superconductor {hujj < 4A). 
Using appropriate materials it should thus be possible 
to reach frequencies from a few GHz to 1 THz, making 
our source compatible with the energy scales of several 
other quantum systems including semiconductor quan¬ 
tum dots. The high frequencies would also allow for sin¬ 
gle photon envelopes on the cm length scale, facilitating 
experiments on quantum non-locality within a cryostat. 

The single-photon source we consider is shown in fig¬ 
ure [TJ It consists of a Josephson junction biased at 
voltage V and embedded in an electromagnetic environ¬ 
ment characterized by the impedance Z{u}) as seen by 
the junction. The impedance is engineered to be high at 
the desired frequency luq, opening a window for photon 
emission at this frequency, as well as high at zero fre¬ 
quency (^(0) > /i/4e^ = Rq). This 0-frequency peak 
is further characterized by a capacitance C, in the sim¬ 
plest case the junction capacitance Gj. When a Cooper- 
pair tunnels across the junction, it gains the energy 2eU 
from the voltage source. However, the electrostatic en¬ 
ergy of the capacitor (G) also increases by the charging 
energy AEc = 2e^ jC. For operation, the released en¬ 
ergy is made to match the frequency of the mode wq, i.e. 
2eV — AEc = fuoQ. After each photon-assisted tunnelling 
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FIG. 1: (a) The single-photon source we consider consists of 
a volt age-biased Josephson junction (JJ) and a parallel junc¬ 
tion capacitor Cj, in series with a semi-infinite transmission 
line (TL) with an impedance step, Zq < R. The TL provides 
a A/2-type standing wave between the step and the junction, 
and results in similar impedance as in (b-c). (b) The desired 
impedance as seen by the JJ, presents a peak at zero frequency 
described by a zero frequency resistance R > h/Ae^ — Rq and 
a capacitance C defining its bandwidth, as well a peak at fi¬ 
nite frequency cuo opening a window for photon emission. A 
photon can be emitted by Cooper-pair (CP) tunnelling when 
the voltage V is chosen so that 2eV = 4e^/2C -I- hujo. Fur¬ 
ther CP tunnelling events are momentarily blocked because 
more energy would be needed to add a second CP to the 
effective capacitor C. After a time RC the capacitor dis¬ 
charges and the next photon-assisted CP tunnelling can oc¬ 
cur. (c) An experimentally feasible realization allowing for 
large zero-frequency impedance R ':$> Rq while maintaining 
compatiblity with standard 50 Q TLs. Here the voltage bias 
is applied to the junction via a, Z = 50 Q TL and a A/4 seg¬ 
ment with Zo > 50 Q. The same TL is used to collect the 
high-frequency response by splitting the signal in high- and 
low-frequency components. In this case the large resistance 
R can be realized on chip as a thin-film resistor. 


event, the high zero-frequency impedance {R > Rq) and 
capacitance lead to a long charging time RC of the junc- 
tion, which temporarily blocks further tunnelling 
[s^. This results in photon anti-bunching in the field 
emitted into the transmission line. In this picture, the 
slow recharging dissipates the rest of the energy supplied 
by the voltage source {2eV — huJo = 4:Ec) in the form of 
a large number of emitted low frequency photons, that 
in an experimental setup are absorbed by the physical 
resistor. 


Methods To quantitatively test if such an interplay 
between the low- and high-frequency parts of the same 
electromagnetic environment is possible, we take a spe¬ 
cific implementation of the environment response func¬ 
tion Z^uj), as described in figure [ija), and study the 
solution of the propagating continuous-mode flux field 
in the neighborhood of the Josephson junction 
{x = 0) and in the transmission line. In the semi-infinite 
transmission line the quantized field can be presented as 
a sum of incoming and outgoing waves (31M^ , 



( 1 ) 




Here, R = sJUJC' is the characteristic impedance, ex¬ 
pressed via inductance L' and capacitance C per unit 
length, and = w\/ CL' is the wave number. The 
incoming (outgoing) wave at frequency w is created by 
the operator a|jj(a;) (dQ^tl^)) annihilated by din(w) 
(aout(w)). They fulfill the standard commutation relation 


ain(w),a|„(a;') 


5{uj — Lo'). A similar solution exists 


also in the cavity region (impedance Zq < R). 

The interaction between the radiation field and 
Cooper-pair tunnelling across the Josephson junction is 
described by the boundary condition (x = 0 corresponds 
to junction location) 


Cjl>(0,t) 


1 d^x,t) 
Lq dx 


=n = L sin 




UJjt - 


( 2 ) 


Here Cj is the junction capacitance, Lq the inductance 
per unit length in the region between the junction and 
the impedance step. The Josephson current is limited 
by the critical current R and controlled by the phase 
(^(t) = 27r<l(0, t)/<I’o. Equation ([2]) corresponds to cur¬ 
rent conservation at the junction. The boundary condi¬ 
tion at the impedance step is linear and can be solved by 
Fourier transformation |35l| . These two conditions can 
now be used to solve the free-space out-field aout(w) as 
a function of the equilibrium in-field ain(w). The general 
solution for the out-field can be written formally as 


dout(^) — 


X 


r{uj)ain{uj) + iIcA{u}] 

V nujir 

U\t,—oo)sm ojjt — ^o{t) 



-oo). (3) 


Here A{uj) (see Supplemental Material) is related to 
the impedance as seen by the Josephson junction, 
Re[Z(a;)] = Zo|A(a;)p, and r{uj) = A{oj)/A{uj)* corre¬ 
sponds to the phase shift in the out-field when Ic = 0. 
The solution is expressed via time-evolution (op erator) 
of the current across the Josephson junction [371 - 1^ . 


17(t,to) = Texp ^ dt' 


Ej cos 




( 4 ) 
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Here Ej = {h/2e)Ic is the Josephson coupling and T 
stands for time-ordering. This is an expansion in terms 
of the phase at the Josephson junction in the absence of 
the tunnelling current, 

Mt) = r + H.c. (5) 

Jo y/uJ 

From here on, we make the natural assumption that 
temperature is low compared to the mode frequency, 
k^T hujo. The flux density of photons due to Cooper- 
pair tunnelling can be evaluated to the leading order in 
the critical current Ic (^1^1^ . 

/■oo 1 

/(w) = J (6) 

= p(2eF-ftw). 

The photon flux is a function of the well-known proba- 
bility density [2^, 



/ OO 1 

(7) 

where J{t) = ([$o{t) — 0o(O)]<Ao(O)^ is a measure of equi¬ 
librium phase fluctuations and is a function of input 
impedance and temperature [^ . The function P{E) de¬ 
scribes the ability of the transmission line to absorb an 
energy E when a Cooper-pair tunnels. In the setup we 
consider, the P{E) function has a simple analytical form 

P{E)^{l-p)PcB{E)+pPcB{E-huJo), ( 8 ) 

where p <C 1 (Supplemental Material) and Pcb.{E) is 
the probability distribution for a high-Ohmic impedance, 
R = Zo^ Rq with cut-off 1/77(7, 

Pcb(E) = ‘ 

(jy/ZTT 

Here, the broadening a = y/SEckuT (assuming IjRC < 
kBT/h) describes thermal fluctuations of the junction 
voltage induced by the series resistor. These fluctuations 
are accounted for as interaction with (thermal) photons 
in the transmission line, treated to all orders. The total 
P(£')-function has not only a peak at E = AEc, but also 
at E = AEc + fiwo, that corresponds to the possibility 
to emit the energy hujo to the lowest cavity mode and to 
use the rest to charge of the junction capacitor. 

Results We now investigate the probability to observe 
two photons with a time separation r. This is usually 
quantified in terms of the second order coherence function 
d^. To calculate this, we need the first order 
coherence function, defined via the photon flux as 


FIG. 2: Second-order coherence based on equa¬ 

tions (llOmill (solid lines) and on semi-analytical approxima¬ 
tion (1121) (dashed lines), at different temperatures. When bi¬ 
ased optimally (the black arrow) we observe anti-bunching 
within RC time-scale and limited by temperature. The 
expected voltage immediately after the first tunnelling is 
pointed by the green arrow. The long-time behavior is re¬ 
produced by the semi-analytical formula, based on calculat¬ 
ing the leading-order emission flux with a slowly evolving 
voltage. At short times, thermal fluctuations increase the 
bunching and wash out anti-bunching approximately when 
T = 60 mK ss 3Ec/kB- (Other parameters are given in the 
Supplemental Material.) 

w ~ wj — 4:Ec / h. The inverse of its bandwidth W deter¬ 
mines the temporal resolution of the photon detection. 
The expression for the unnormalized second-order corre¬ 
lation function of photon detection reads 

G(2)(r)= (9) 

using the shorthand notation = H^ duJiF(LJi) 

and = aout(wi). The dominating contribution comes 
from the fourth order in R (Supplemental Material), 

^times ^ '' '' '' ' 


Dfc 7*00 

= — I dujP-^u;\E{io)\^f{oj). 

Here, E(uj) describes the measurement band-pass fil¬ 
ter [4l|| . centered at the chosen measurement frequency. 


where = A*{ijJ 2 )A{wo)A(w 4 ) 6 (ijJi -b W 2 — wa — 
a; 4 ) and = H^ J^^dp. The operator = 


giivt gjj-p 


i[$oit) - LOjt] 


can be interpreted as a creation 


of a photon with frequency w via Cooper-pair tunnelling 
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at time t. The final expression for the second order co¬ 
herence can be presented as 




Gfiir) 


= 1-b 


g2nd(^) 

<^d(0) 


+ G{r), (11) 


where G{t) captures correlations between consecutive 
Cooper-pair tunnelling events and the other terms de¬ 
scribe the emission from uncorrelated Cooper-pair tun¬ 
nelling. When the temporal resolution of the detection 
cannot resolve individual tunnel events, i.e. for W —>■ 0, 
we get t/(r) —>• 0 and 2, which is characteristic 

for chaotic (thermal) light, having no phase coherence 
between different photons [i^. For times much longer 
than the temporal resolution r ^ W~^ and recovery time 
T ^ RC, only the first term is finite leaving (r) = 1, a 
characteristic of uncorrelated (Poissonian) photon emis¬ 
sion. Outside these limits, we see clear effects of G{t), as 
shown in Fig. [51 where we plot the second-order coherence 
calculated numerically using equations (IIOIIIII) . We con¬ 
sider a cavity length such that a resonance mode appears 
at aJo/27r = 5 GHz, for R = 4i?Q and ZojR = 1/10. The 
detection filter is centered at the mode frequency with a 
bandwidth of 1 GHz. The P(i?)-function (inset) is close 
to the analytical form in Eq. (|5]) and the simple interpre¬ 
tation presented in figure [I] implies clear anti-bunching 
[g(^)(0) < 1] of the outgoing radiation. Indeed, when bi¬ 
ased optimally at 2eV = AEc + tEjo (black arrow in the 
inset), we observe a clear anti-bunching over a time-scale 
given by the RC-time. The depth of the anti-bunching 
dip is limited by temperature, and increasing it beyond 
the charging energy, /cbT > SEc, the output radiation 
instead becomes bunched > !]■ 

Classically, at zero temperature, a sudden 2e-reduction 
of the equilibrium capacitor charge (Q = CV) induces 
recharging dynamics v{t) = V — (2e/C')e“*/^‘". A first 
guess would be that (in the full quantum treatment) the 
probability to detect a second photon after a time r is 
given by the photon emission rate at voltage v(t). Fol¬ 
lowing this simple idea, we write down a semi-analytical 
formula for the expected second order coherence. 


( 2 )( ^ Ibw [^(^)] P [2eu(r) - feu] 

^ /3^Re[Z(w)]P(2eF-M ’ ^ ^ 

In figurelU we see that equation ([T2]) reproduces the long¬ 
time behavior, whereas can fail when the separation be¬ 
tween the tunnelling events becomes short, i.e. when the 
tunnelling processes start to overlap. This can also be 
understood more mathematically, as detailed in the Sup¬ 
plemental Material. The form implies that to suppress 
g^^^{0) at the optimal bias point, we need the width of 
the P{E) to be smaller than the bias jump iEc, which 
leads to the conditions k^T Ec and R ^ Rq. 

For further illustration of recharging effects in this sys¬ 
tem, we study the dependence of bunching on the bias 
voltage, see figure [3] For high bias voltages the second- 
order coherence shows super-bunching (g^^^(O) > 2), 




Frequency (GHz) 


FIG. 3: (a) Second-order coherence for different bias 

voltages, based on equations (llUmil) (solid lines) and on ap¬ 
proximation (1121) (dashed lines), for the system as in figure [2] 
with T = 10 mK. (b) The P(i?)-function and the positions of 
the studied voltage-bias points, marked by the arrows of the 
corresponding colors, (c) The photon-flux density as a func¬ 
tion of the bias voltage in the neighborhood of the studied 
mode at a;o/27r = 5 GHz. Depending on the bias point, the 
second-order coherence can have a deep antibunching dip, a 
local maximum, or show a steep decay, as a function of the 
measurement separation-time r. The steep decay or local 
maximum occurs since above the optimal bias (black arrow) 
the first emission rate is decreased, while the rate for the sec¬ 
ondary emission can be always higher, or sweep through the 
maximum, during the recharging. 


whereas for lower bias voltages both anti-bunching and 
bunching can co-exists for a single curve. These prop¬ 
erties are qualitatively reproduced by equation (HU, and 
occurs here since the first emission rate is decreased com¬ 
pared to the optimal point {2eV = AEc + Swq), while 
the rate for the secondary emission is always higher, 
or sweeps through the maximum, during the recharg¬ 
ing. We also observe that g^^^(O) can be suppressed 
below the value of the optimal point, by a bias-voltage 
2eV < AEc + kujQ, with the cost of a reduced photon flux. 
This is a result of a reduction in thermal-fluctuation trig¬ 
gered two-photon emission processes. 

In conclusion, by extending the well-known P{E)- 
theory to two-Gooper-pair processes (i.e. to 4th order 
in Ej) we have shown that the Goulomb blockade of 
Gooper-pair tunneling allows for creation of strongly 
anti-bunched microwave photons from a simple dc bias 
under experimentally realistic conditions. The validity 
of our theoretical approach is limited to tunnelling rates 
lower than IjRC, but we expect anti-bunching to sur- 
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vive also at higher tunnelling rates. This regime would 
be characterized by modified junction voltage dynamics, 
and is indeed an interesting regime to explore, both ex¬ 
perimentally and theoretically. The proposed source of 
antibunched photons can be transformed into a photon 
gun by replacing the Josephson junction by a dc-SQUID. 
Cooper-pair tunnelling can then be efficiently suppressed 
by threading a half flux quantum through the SQUID and 
a flux pulse short compared to RC will then release a sin¬ 
gle photon on demand with high probability for large Ej. 
Therefore, this system allows for an on-demand photon 
emission at a very high rate and in a wide range of fre¬ 


quencies. Such a bright microwave single photon source 
would be useful in various microwave quantum measure¬ 
ment setups because of its perfect amplitude squeezing 
and well defined power. 
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Supplemental Material 


Factor A{lo) 


Fourier transformation of the boundary conditions at the Josephson junction and at the impedance step leads to 


A{u}) 


O /-^o ^ — 2ik^ d 


(1 + ^) + (^ 


l)C*(a;)' 


(SI) 


Here C(a;) = l — iZ^Cju}, is the wave number, and d is the length of the cavity. A detailed mathematical 

derivation of A{uj) is given in reference [35] (where the corresponding function is called A(a;)). 


P(i7)-function 


The real part of the impedance seen by the tunnel junction can be deduced to be Re[2'(u;)] = Zo|A(w)p (see 
reference [35]). In the case Zq ^ R the transmission line has (A/4-type) resonance frequencies approximately at 
(jJn = (2n -I- l)u;r, where Wr = 7r/2dY^C'oTo and n £ [0,1,2,...]. Near each resonance the real part of the tunneling 
impedance is approximately a Lorentzian. In the neighbourhood of the lowest mode wq we can write (keeping the 
notation similar with reference [25]), 


Re[Z(w)] 


1 r 

a 1 + 4(W - Wo)2r2 


2 Q 


■S{uj - Wo). 


(S2) 


Here Ct = 7r/4wo.^o a-nd F = ttZq/ARluq. The corresponding P(i?)-function, in the d-function approximation and at 
zero temperature, has the form [25] 


OO ^ 

P{E) = e-P V ^6(E - khujo). (S3) 

k\ 

Here p = l2Ct)fhujQ compares the effective capacitive energy of a single Cooper pair to the mode frequency. In 

the case Zq ^ R we have p = AZq/Rq. For a finite linewidth of the mode the resulting P(i5)-function broadens 
accordingly. 

In the case considered in this article, Zq <C R-, the transmission line has (A/2-type) resonance frequencies approxi¬ 
mately at nwr, where Wr = "k j and n £ [0,1,2,...]. We have therefore a peak also at zero frequency. We can 
take here the approximation C(w) ~ 1 which leads to 

Re[Z(a;)] « - ^ (S4) 

1 -I- 

Zo ujo / 

This means that the cut-off frequency of the zero-frequency peak is IjRC, where C = ttIZquiq. For the used parameter 
range this effective capacitance dominates the junction capacitance (C « 2C'j). 

The effect of the junction capacitance Cj is more pronounced for the first resonance where it induces a shift to a 
lower frequency and a reduction of the impedance. We find an approximative correction to the resonance frequency 
in the considered parameter region. 
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(S5) 


Here we use the notation Wc = 1/RCj. In the considered case we have wo/wr ~ 5/7. This allows for an analytical 
expression for the impedance in the neighbourhood of wq. However, the form of the Lorentzian in this case is rather 
cumbersome and less useful. According to our numerical calculations, the corresponding probability parameter satisfies 
p <C 1. This can also be deduced from the relative heights of the first two peaks in the P(A)-functions plotted in 
figures 2-3, bearing in mind equation (S3). This observation, and that the total P(P)-function of a sum of two tunnel 
impedances is a convolution of the two P(£')-functions obtained for individual impedances, leads to equation (8). 
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Second-order coherence function 

Two-photon emission is characterized by the second-order coherence, 

G(2)(r)= . 

An expansion to the leading-order (second-order) in the junction’s critical current describes photon emission from 
single Cooper-pair tunneling. The straightforward calculation is presented in reference [35], and involves the summa¬ 
tion. 


(a’*’a’''aa) —!► ^ , (S6) 

done with the constriction i + j + k + I = 2. The sub-indices correspond to the order (in the critical current) of 
operators a. When z = 0 oiW = 0, we have a contribution proportional to the Bose factor at the measurement 
frequency, which can be neglected in the frequency range we are interested in, Huj ^ feeT. The main contribution 
comes when i = I = 1 (and j = fc = 0), giving the unnormalized second-order coherence, 

P[fi(tcj - - cc')]. 

Jbw 

This is characterized by the probability P(hujj — Huj — Huj'). This means that in the leading order the observed photons 
originate from single CP tunneling processes, each of them releasing an energy hujj, and that this energy has to be 
split into two individual photon energies. At the frequency range we are interested in, ui ~ wq, such photon pair 
production is negligible as long as k^QT/h <C cug- Therefore, it is the fourth-order contribution in the critical current 
[equation (10)] that is dominating at low temperatures. A derivation of this contribution is given in reference [39]. 


Semi-analytical formula 

For an analytical expression of the second-order coherence, we consider the limit of a wide bandwidth (high temporal 
resolution), and that the field operators 


a(r)- [ e*“(‘+")A(o.)A,^, (S7) 

JBW 

can be treated local in time (discussed in more detail below). This means that in the following ti < t 2 and 

We then study contraction of outer and inner tunneling operators in equation (10). We use the property of bosonic 
pair correlations, 


n, 


exp 






'j = exp 

-'^ninjJ{ti - tj) 
i<j 


where rii is an integer 1 or -1, and = 0- We get, 


Gith('^) oc 


(< 2 -I-t) giiiia (ta-l-T) ^21^4*4 


«/BW times 

^ ^ t4) + ^(t2—t3)g^PI 


(S8) 


Here we have introduced a function that describes interaction between the contraction pairs, 


Jpi = J[ti - ts) -I- J(t2 - t4) - J{ti - 12) - J{h - G). (S9) 

If Jpi = 0, the pairs decouple and the result of the time integration is proportional to a product of two P(£’)-functions, 
P(hujj — huj) X P(hujj — Huj'). Here uj = uji = UJ 4 and uj' = UJ2 = ^3- This decoupling occurs always for long r (which 
means here long <2 ~ ti or equivalently long ts — t^). It can also occur for intermediate r with a modified energy 
argument, as demonstrated below. 
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Let us consider the case Jpi 7 ^ 0 and when the most relevant contribution to the second-order coherence comes 
from short-time behaviour inside the pairs, |t| = \ti — ^ 4 ] < r and \t'\ = \t 2 — < r, i.e. there is a fast decay of the 

term increasing t or f. This helps us because the time-ordering between the pairs is fixed 

and the integrations between the pairs do not overlap (as assumed at the beginning, ti < t 2 and > ti). When 
k-BT/h < \/RC, the time-scale of the decay is of the order, or less than, RC. This means that we are safe (at low 
temperatures) if r > RC. 

In the next step, when estimating the term Jpi, we notice that all the time arguments of the individual phase 
correlation functions are of the order t. This allows us to consider an approximation of the phase-correlation function 
(when evaluating Jpi), 


Jit) = c-D'\t\-iV't, (SIO) 

where D\ V', and c are constants for a fixed r. This is a linear approximation for times in the neighbourhood of r. 
The time-dependent imaginary part reads V' { t ) oc V { t ) = (2e/C)(l — corresponding to a momentary value 

during a classical RC-recovery of the voltage after a Cooper pair tunneling event. Inserting approximation (SIO) to 
equation (S9) leads to Jpi = —2iV'it2 — ta), i.e. Jpi is purely imaginary. This term adds to the imaginary part of 
Jit 2 — ts) in equation (S 8 ) and describes a shift in the junction voltage V ^ V — V. Integration over the times and 
frequencies then lead to approximation ( 12 ), 

_ /bw [■^(^)] P ('^) ~ ('Slli 

^ /3^Re[Z(a.)]P(2eR-M ’ ^ ^ 

In the case of a high temperature, k^T/h ^ 1/RC, the phase-correlation function J(t) at intermediate times is 
quadratic rather than linear. This leads to the failure of the formula in a wider range of time r as for k^T jh ^ Ij RC. 
The general failure at short times can also be understood as a neglection of correlated voltage fluctuations between the 
Cooper-pair tunnelings, and as a neglection of two-photon emission via 4e-tunnelling (’’virtual” initial Cooper-pair 
tunnelling). The formula also neglects interference effects between photons that lack relative phase memory, which 
increases 5(2)(0). 


Parameters in the numerical simulations 


In figures 2 and 3 we use a free-space transmission line R = 4i?Q and an impedance step Zq/R = 1/10. Here 
Rq = h/4e^ ~ 6.5 kfl is the superconducting resistance quantum. The resonance frequency (a;Q/27r) is at 5 GHz 
and the junction capacitance is Cj = 25 fF. The measurement band-pass filter P(a;) was chosen to be a Gaussian 
with a full width at half maximum 1 GHz. As discussed above, the effective capacitance C differs from the junction 
capacitance, and becomes C = 50 fF. The value is determined from the RG recovery of the imaginary part of J{t). We 
have therefore a charging energy Eq = 1-6 ^eV (fes x 18.6 mK) and a charging time RC = 1.3 ns. Equations (lO)-(ll) 
were evaluated numerically using a three dimensional fast Fourier transformation. Due to a finite discretization of 
the time and frequency spaces (2^^ points, where N = 9), the plotted values of can have a small numerical 

error, which we estimate to be limited by ± 0 . 02 . 


Feasibility of the voltage bias and high resistance environment 

Regarding the experimental feasibility, one question is the level of control of the bias voltage and temperature. 
Thermal fluctuations in the theory of figure 1(a) are accounted for by the model as it treats the interaction with 
bosonic modes in the bias line to all orders. The experimentally relevant setup, figure 1(c), behaves in the same way, 
but here the effective temperature has two contributions; one from the on-chip resistor R, and one from fluctuations 
in the applied voltage. Using low-impedance cold voltage dividers, the fluctuations of the applied bias voltage can be 
made as low as a few 100 nV (see reference [28]), which is small compared to the voltage fluctuations of the on-chip 
resistor at 10 mK (several gV). Therefore the relevant difficulty is the cooling of the on-chip resistor. 

The on-chip resistance R > Rq can be realized, for example, using a thin-film chromium wire. Thermalisation of 
on-chip resistors is rather well established and is discussed in detail in Phys. Rev. B 76, 165426 (2007). We explore in 
figure 2 how cold the resistor must be for our proposal to work and obtain experimentally achievable temperatures. 
The resistor will have stray capacitances, but in the acceptable range for the proposal to work. 
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Feasibility of the on-demand photon emission 

In the last part of the article we describe how one can realize a photon gun by using a SQUID rather than one 
Josephson junction. By changing the magnetic flux through the SQUID one can reach the right regime of photon 
emission for a controlled amount of time. One concern here could be that hysteresis in the junction dynamics would 
hinder the relaxation to the desired voltage bias condition. Hysteresis is avoided since the flux pulse is used to tune Uj 
to zero, so that there will be no nonlinear effect due to the junction which could lead to hysteresis. If this flux value is 
maintained for a time long compared to the relaxation-time of the electromagnetic environment, that is here i?C, any 
sub-gap voltage can be reached. The increase (and decrease) of Uj should be fast compared to RC but slow compared 
to l/wo, not to create photons by the pulse itself. If Uj is increased slowly compared to RC instead, the junction 
voltage can latch to 0 after a photon has been emitted, due to the possibility of consecutive Cooper-pair tunnelings 
with no photon emission. This effect can actually be beneficial for our proposal, as it deepens the anti-bunching 
(but goes beyond the theory developed in this paper). We also note that sub-gap quasiparticle tunneling will not be 
induced, when the applied flux-pulse is slow in the timescale of the (inverse) superconducting gap. 


